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THE RANGE OF MULTIPLICATIVE FUNCTIONS ON C[x],M[x] AND Z[x] 



CHRISTOPHER D. SINCLAIR 



Abstract 

Mahler's measure is generalized to create the class of multiplicative distance functions. These functions measure 
the complexity of polynomials based on the location of their zeros in the complex plane. Following work of 
S.-J. Chern and J. Vaaler in pQ, we associate to each multiplicative distance function two families of analytic 
functions which encode information about its range on C[x] and WL[x]. These moment functions are Mellin 
transforms of distribution functions associated to the multiplicative distance function and demonstrate a great 
deal of arithmetic structure. For instance, we show that the moment function associated to Mahler's measure 
restricted to real reciprocal polynomials of degree 2iV has an analytic continuation to rational functions with 
rational coefficients, simple poles at integers between — N and N, and a zero of multiplicity 2N at the origin. 
This discovery leads to asymptotic estimates for the number of reciprocal integer polynomials of fixed degree 
with Mahler measure less than T as T — > oo . To explain the structure of this moment functions we show that 
the real moment functions of a multiplicative distance function can be written as Pfaffians of antisymmetric 
matrices formed from a skew-symmetric bilinear form associated to the multiplicative distance function. 



1 . Introduction 

This manuscript is concerned with measures of complexity of polynomials which respect 
both the algebraic structure and topology of C[x] (as generated by all open sets in all finite 
dimensional subspaces of C[x]). As such, we are interested in functions from C[x] to the 
non-negative reals which are continuous (as a function on coefficient vectors) on all finite 
dimensional subspaces of C[x] and behave nicely with respect to multiplication and scalar mul- 
tiplication. The most important requirement of these functions is that they be multiplicative. 
As we shall see, multiplicativity is a very strong condition which allows for many interesting 
theorems. 

The following axiomatization suggests itself: A function <I> : C[x] — -> [0, oo) will be called a 
multiplicative distance function if 

Al. $ is continuous, 
and for all f,g € C[x] and aieC, 

A2. $ is positive definite: <&(/) = if and only if / is identically zero, 

A3. <I> is absolutely homogeneous: §>(wf) — \w\$>(f), and 

A4. <f> is multiplicative: &(fg) = $(/)$(.g). 
The nomenclature stems from the fact that multiplicative distance functions restricted to finite 
dimensional subspaces of C[x] are distance functions in the sense of the geometry of numbers. 
We will refer to <&(/) as the ^-distance of / to the origin or simply the distance of /. 

It is easily seen that $ is uniquely determined by its action on monic linear polynomials. 
That is, if 

N 

f(x) =al[(x~ ln ) : (1.1) 
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then there exists a function : C — > (0, oo) such that 

N 

* (/) = H I] 0(7n). 

n=l 

The function will be known as the root function of <i>. The best known example of a 
multiplicative distance function is Mahler's measure, denoted by p, and defined by 

N 

Kf) = m n max { 1 ' i^i}- 

From the definition of p it is clear that Mahler's measure satisfies Axioms A2, A3 and A4. It 
is less clear that Mahler's measure is continuous on C[x], but this was proved by K. Mahler in 
1961 [5]. 

Mahler's measure can be extended to the algebra of Laurent polynomials by using 1 = p(l) = 
/z(xx _1 ) to write p{x~ r ) = = 1. If <5>(x™) = 1 for every n > (or what amounts to 

the same thing, 0(0) = 1) then we shall say $ is shift invariant. Shift invariant multiplicative 
distance functions can be naturally extended to the algebra of Laurent polynomials by setting 
(pix- 1 ) = 1. 

The continuity of $ controls the asymptotic behavior of cf>. And in fact, this asymptotic 
condition produces a classification of multiplicative distance functions. 

Theorem 1.1. Suppose that $ is a multiplicative distance function with root function <f>. 
Then, 

<t>{i) ~ hi as M -> °°. 

Conversely, if if) : C — > (0,oo) is a continuous function such that t/>(7) ~ ItI as ItI ~^ °°> then 
ip is the root function of a multiplicative distance function. 

1.1. Examples of Multiplicative Distance Functions 

Theorem 11.11 gives us a way of producing examples of multiplicative distance functions, 
and in this section we will introduce another method for constructing multiplicative distance 
functions. 

A Laurent polynomial g{x) g C[x,l/x] is said to be reciprocal if g(l/x) — g{x), and 
the algebra of reciprocal Laurent polynomials is given by C[x + 1/x]. Clearly the algebra 
of reciprocal Laurent polynomials is a subalgebra of C[x, 1/x] and hence we may speak of the 
Mahler measure of a reciprocal Laurent polynomial. We define the reciprocal Mahler's measure, 
p, of / € C[x] to be the Mahler's measure of the reciprocal Laurent polynomial f(x + 1/x). 
That is, p(f) = p{f{x + 1/x)). It follows that the root function of p is given by 



7 i — ► (i(x + 1/x — 7) = max |l, 

This definition is independent of the branch of the square root used and is easily seen to satisfy 
the conditions of Theorem 11.11 

The procedure used to create the reciprocal Mahler's measure may be repeated to create 
multiplicative distance functions formed from Mahler's measure restricted to other subalgebras 
of C[x, 1/x]. In particular, if G(x) G C[x, 1/x] is a fixed Laurent polynomial we may create 
a multiplicative distance function by considering Mahler's measure restricted to C[G(x)] C 
C[x, 1/x]. Thus we define G* p : C[x] -> [0,oo) by G*p{f) = p(f o G). It can be verified that 
G* p satisfies all the axioms of a multiplicative distance function. The notation for G* p stems 
from the fact that if we view G as the natural map C[x] — > C[G(x)] then G*p is the pullback 



max < 1 . 
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C[x] — C[G(x)} 
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Given < t < 1 we define the t-reciprocal Mahler's measure, m, to be the pullback of \i 
through the Laurent polynomial x + t/x. In this context [i\ is the reciprocal Mahler's measure 
and /iq is Mahler's measure. Thus, as t varies from to 1 we have a 'path' of multiplicative 
distance functions whose end points are fi and p. 



1.2. Potentials and Jensen's Formula 

The prototypical multiplicative distance function, Mahler's measure, satisfies an important 
integral identity. If / is given as in (jl.ip then Jensen's formula implies that 



/V 



/"(/) = M Y[ m ax{l, 



|7„|} = exp<{ — 



2tt 



log|/(e je )| d6 



(1.2) 



The right hand side of this equation is an example of an equilibrium potential. By generalizing 
the right hand side of (|1.2p we may produce examples of multiplicative distance functions 
which are associated to compact subsets of C. Multiplicative distance functions of this sort 
were considered from the standpoint of equidistribution by R. Rumely in [7]. 

Let if be a compact subset of C and let v be a probability measure whose support is contained 
in K . The potential of v is defined to be the function p v : C — > [0, oo) specified by 



Pu(l) = exp 



log \z — j\ dv(z) 



K 



It is a fundamental result of potential theory that p v is upper semicontinuous. Moreover p v (7) ~ 
I7I as I7I — > 00 and thus if p v is in fact continuous then it is the root function of a multiplicative 
distance function. 

If we denote the set of probability measures whose support lies in K by M{K), then under 
fairly mild conditions on K there is a unique probability measure vk £ M(K) which minimizes 



1(1/) 



log[p„(7)| dv(p() over all v € M(K). 



K 



For instance, the minimizing measure is unique if there exists at least one v € M(K) with 
< 00. When exists this measure is known as the equilibrium measure of K and the 



quantity c{K) = e 



is known as the capacity of K. We will denote the potential of vk 



simply by pk ■ This potential is called the equilibrium potential of K . 

If K is regular with respect to the Dirichlet problem then K has positive capacity and px is 
continuous. In this situation px is the root function of a multiplicative distance function which 
will be denoted Pk- For instance, if if is a simply connected compact subset of C which does 
not consist of a single point then K is regular with respect to the Dirichlet problem and we 
may speak about the multiplicative distance function Pk- Explicitly, 



N 



PkU) 



\\PK{ln) 



exp 



log \f(z)\ dv K (z) 



For example, Mahler's measure can be represented as Pd where D is the closed unit disk. 
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It is a well known fact of potential theory that pk{i) > C {K) with equality if and only if 
7 £ K. Of particular importance are multiplicative distance functions associated to simply 
connected compact sets of capacity 1. In this situation if K contains the origin then pk(0) = 1 
and hence Pk is shift invariant. 

As the next theorem demonstrates, there is a strong connection between multiplicative 
distance functions formed from certain compact sets K and those formed by the pullback 
of Mahler's measure by certain rational functions. 

Theorem 1.2. Let q(x) £ C[x] be a monic polynomials of degree M and define G(x) = 
q(x)/x M ~ 1 . If G(x) is a conformal map from C \ D onto its image, then G*/i = Pk where K 
is the complement in C of G(C \ D). 

Theorem ll.2l is well-known to experts in potential theory (though perhaps not in the language 
used in this manuscript). From the definitions of G*/i and Pk, the equation G*fi = Pk may 
be thought of as an analog of Jensen's formula. 

Corollary 1.3. Let < t < 1 and define E t C C to be the compact set given by 

( x 2 y 2 ) 

and define 

E 1 = {x + ty:y = Q,x£ [-2,2]}. 
Then, for any t £ [0, 1], the t-reciprocal Mahler's measure, fi t , J 's equai to PE t - 

As t varies from to 1, E t deforms from the unit disk through a series of regions bounded 
by ellipses to the degenerate ellipse given by the interval [—2, 2] on the real axis. All of these 
compact sets have capacity 1. Thus our 'path' of multiplicative distance functions formed from 
the pullback of /i through x + t/x as t varies from to 1 can also be thought of as a 'path' 
of shift invariant multiplicative distance functions formed from the family of ellipses E t as t 
ranges over the same values. 

1.3. Star Bodies and Distribution Functions 

By identifying each polynomial of degree N with its vector of coefficients, the set of poly- 
nomials in C[x] of degree N may be identified with the vector space C N+1 . To each a £ C N+1 
we define the polynomial a(x) by 

N+l 

a(aO = J2 *nX N+1 - n . 

n=l 

We may regard $ as a function on C N+1 by setting ^(a) = $(a(a;)). As such $ satisfies all the 
axioms of a vector norm except the triangle inequality. The 'unit ball' of $ is thus not convex. 
That is, the set 

Vjv($) = {ae C N+1 : $(a) < 1} 

is a symmetric star body about the origin which will be referred to as the degree iV complex 
unit star body of Similarly the degree N real unit star body is defined to be the set 

U N ($) ={ae R N+1 : *(a) < 1}. 

The absolute homogeneity of $ implies that the set of polynomials of degree N in C[x] with 
distance bounded by T > is the dilated star body TVn- 
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As a first application of the theory of multiplicative distance functions, S-J. Chern and 
J. Vaaler devised a procedure for determining the volume (Lebesgue measure) of Wj\r(/Lt) and 
then used this to give the main term in an asymptotic estimate for the number of polynomials 
in Z[x] with degree at most N and Mahler measure bounded by T as T — > oo pQ. Their idea is 
more generally valid, and we will give similar estimates for the reciprocal Mahler's measure. 

Theorem 1.4. Let $ be a multiplicative distance functions. Then, as T — > oo, 
#{ae Z N+1 : *(a) <T}= vol(U N ($))T N+1 + 0(T N ). 

Proof. See H Ch. VI, §2] or [TJ §12] □ 

In order to determine the volumes of £//v(/i) and Vjv (/•*), Chern and Vaaler introduced 
two families of analytic functions which encode information about the range of values of u. 
restricted to polynomials with real and complex coefficients. Their techniques generalize to 
other multiplicative distance functions and the analogous analytic functions demonstrate a 
great deal of structure which can be used to learn information about the range of values of a 
multiplicative distance function. 

We define the degree N monic restriction of $ to be the function $ : — > (0, oo) given by 

$(b) = $ + b n x N - n ^ . 

That is, $ is simply $ restricted to the set of (non-leading) coefficient vectors of monic 
polynomials of degree N. We use Xn and X2N to denote Lebesgue measure on Borel subsets of 
M. N and (respectively) and define the distribution functions /iv,/ijv : [0, 00) — > [0, 00) by 



/ jv ($;C) = A A r{beM Ar :$(b)<e} 



and 



1 N 



(<J>;e) = A 2A r{beC Ar :$(b)<e} 



By identifying R w with the set of monic coefficient vectors in R JV+1 , /jv($;£) is simply the 
volume of the intersection of the dilated star body £Wjv($) with M. N . In this way /at encodes 
information about the range of values $ takes on monic polynomials of degree N in R[x]. 
Similarly /ijv($; £) encodes information about the range of values $ takes on monic polynomials 
of degree N in C[x]. For instance, the volumes of Un{^) and VV(<E>) can be discovered from 
and M$;0- 

Theorem 1.5. The supports of /jv($;£) and /iat(<E>;£) are bounded away from 0, and as 
£ -> 00, 

- 0(( N ) and h N {fr,0 = 0(£ 2Ar ). 

Moreover, 

lim /jV ^ ;C) = Ajv(fcfo-i(*)) and lim fejv( * ; = \ 2N (V N -i(3>)). 



The Mellin transform of these functions is then given by 



? 



? 



where s is a complex variable. From the asymptotic formulae for /jv($;£) and /i_/v(<I>; £) it is 
easy to establish that the integral defining /jv($; s) converges when 5R(s) > TV, and the integral 
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defining /ijv($; s) converges when 3?(s) > 2N. Moreover, by Morera's Theorem /jv and hisr are 
analytic functions in their respective domains of convergence. These analytic functions encode 
information about the range of values <E> takes on monic polynomials of degree N in R[x] and 
C[x] respectively. For instance, the volume of Un{$) also appears as a special value of /at($; s) 
and similarly the volume of Vjv(<I>) appears as a special value of Hn{^\ s). 



Theorem 1.6. The volume ofU N (<&) is given by 

\ N+1 (U N ($)) = 2fa($;N+l), 
and the volume of Vn($) is given by 

^N + 2(V N m = 2ttM$; 2N + 2). 



Beyond the computation of the volumes of Un(&) and Vat($), any analytic continuation of 
/at($; s) and ft.Ar(<I>; s) beyond the range of convergence may yield further information about 
the range of values of $ which may not be realizable from other methods. 

It should be remarked that the proofs of Theorem 11.51 and Theorem 11.61 do not rely on the 
multiplicativity of <i>. 

It is not obvious that, for any choice of <E>, the integrals defining /at($; s) and hjs[(&', s) can 
be expressed in terms of well-known analytic functions. As a first step in this direction, we 
view the integral defining /jy as a Lebesgue-Stieltges integral and use integration by parts to 
write 



Ms) 



+ 



1 



C s df N (0- 



(1.3) 



It follows from Theorem 11.51 that /n{0) — and that /iv(£) is dominated by C£ for some 
constant C. Consequently, the first term in (|1.3|) is 0. From the definition of d/W(£) we can 
write 



7jv(*;s) = -F N ($;s) where F N {<£;s) 



Similarly, 



h N ($;2s) = — H N (<f>;s) where H N {<f>;s) 
is 



$(b)- s d\ N (b). 



$(b)-^ s d\ 2N (h). 



Both Fjy(&; s) and Hn(&; s) converge to analytic functions in the region 9?(s) > N. We will 
call these the real and complex degree N moment functions of <& (respectively) . 



1.4. Examples of Moment Functions 

Chern and Vaaler's original motivation for computing the moment functions for Mahler's 
measure was provided by Theorem 11.41 Amazingly, their computation revealed that Hn (p; s) 
and Fn(p; s) analytically continued to rational functions of s with poles at integers and a high 
order zero at s = 0. Moreover, they showed that both F/v(/i; s) and tt~ n Hn(/J.; s) have rational 
coefficients. 

Theorem 1.7 (S.-J. Chcrn, J. Vaaler). Fn(h;s) and iJjv(/i;s) analytically continue to 
rational functions of s. In particular, 

N N 

AM ±J ~ s — n 

n—l 
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If M is the integer part of (N - 1) /2 then 

M s M ( 2m \ N ~ 2m 

M«s)=C N l[ s where ^ = 2^(3^+1) M 

m=0 v ; m=l v ' 

This surprising result provides additional motivation for determining the moment functions 
of other multiplicative distance functions. And in fact, the author's original motivation for 
introducing multiplicative distance functions and their moment functions was to create a 
context in which the surprising rational functions identities of Chern and Vaaler could be 
explained. The next result shows that much of the structure evident in F 'n (//; s) and Hn (//; s) 
carries over to the moment functions of the reciprocal Mahler's measure. 

Theorem 1.8. Fn(p; s) and Hn(p; s) analytically continue to rational functions of s. In 
particular, 

%»)=2vn^. (i.6) 

n=l 

If J is the integer part of (N — l)/2 then 

F N ( P ;s) = v N ]l s2 _ {N _ 2jr where v N = -^—) . (1.7) 

A variation of (|1.6p was established in [9]. 

The parity (evenness/oddness) of H N {p;s) and F N (p;s) should be mentioned. This sym- 
metry seems to arise from the fact that p is the pullback of Mahler's measure through the 
polynomial x + l/x. The Mellin transform translates the symmetry x i— > 1/x to the observed 
parity in the moment functions. We may view the parity of Hn(p; s) and Fn(p; s) as a kind 
of functional equation, and it seems likely that the mechanism which produces this functional 
equation will produce functional equations for moment functions for other multiplicative dis- 
tance functions formed from the pullback of Mahler's measure through other rational functions. 

As p and p are the 'endpoints' of a 'path' of multiplicative distance functions so are Fn(ii; s) 
and Fpj(p; s) the 'endpoints' of a 'path' of moment functions, and similarly for Hn(/j,; s) and 
Hjy(p;s). Much of the structure present in Theorem 11.71 and Theorem 11.81 carries over to the 
moment functions of p t f° r < t < 1. By investigating the qualitative properties of the moment 
functions of p t we- may hope to learn how the structure of moment functions relates to the 
underlying multiplicative distance functions, in particular for those moment functions which 
arise as pullbacks of Mahler's measure through rational functions. 

Theorem 1.9. Let < t < 1. Then, F/v(//t;s) and ffjv(/xt;s) analytically continue to 
rational functions of s. In particular, 

n N s N » (l-t 2 »).s + (l + t 2 > 
Hn(pus) = 11 , 

71=1 

and Fpf(p,t', s) £ Q[t](s). Moreover if J is the integer part (N — l)/2 then Fwipt', s) has simple 
poles at ±N, ±(N — 2), ... , ±(JV — 2 J), a zero of multiplicity J at s = and J other real zeros 
on the negative real axis. 



We will leave this theorem unproved since its proof relies on the same methods we will use 
to establish Theorem ll.81 Notice that when t — and t — 1 the formula for Hn{pu s) coincides 
with the formula for Hn(p; s) and Hn(p; s) respectively. It should be remarked that a closed 
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form for Fn(im] s) can be discovered using the same method of proof as Theorem 1 1 . 81 and this 
expression agrees with those for Fn(h;s) and Fx(p;s) when t = and t = 1. However, the 
closed form for _F/v(/ii t ; s) is more complicated than those given for F^(p; s) and F^(p; s) and 
in its place we present Figure [T] 




Figure 1. The location of the zeros and poles of Hq^; s) and Fs(fi t ; s) 



In both cases, as t — » the nontrivial zeros (those not located at s = 0) move to cancel the 
poles located at negative integers. And, since fit is the pullback of Mahler's measure through 
x + t/x, this seems to suggest that the poles at negative integers in H^{pt] s) and Fjs[{fit] s) 
arise from the t/x factor in x + t/x. Note that as t — *■ 1 the nontrivial zeros of Hj^(p t ]s) 
approach ™oo while the nontrivial zeros of F^^fit] s) approach 0. To address this disparity, let 

^ f ) = {bel A ': / I t (b) = l} and V^(^) - {b e C N : ft(b) = 1}. 

The definitions of Hivifit', s) and Fn{pu s) imply that 

Mm F N (fi t ;s) = \ N (U N (fi t )) and lim H N {fi t \ s) = \2N(V N {pt))- (1.8) 

s — >oc s — >oo 

From the definition of fit we see that b is in Wjv(A*t) exactly when + 53 n =i b n x N ~ n has 
all of its roots in the elliptical region As £ — > 1, E t approaches the interval [—2,2] on the 
real axis and the volume of Vjv (/•**) approaches since if b S Vn{p) then in fact b 6 Un{p)- 
If we momentarily identify with R 2Ar then we see that Vjv(p) is a subset of codimension 
N in M. 2N . It is exactly this fact which explains why Hn{p\ s) has 2N poles and only N zeros. 
Moreover the fact that as t — > 1 the non-trivial zeros of Hpf(fj, t ;s) tend toward — oo verifies 
our intuition that the volume of VN(fJ-t) is tends toward 0. On the other hand, since Un(p) has 
positive Ajy-measure we expect F^(p;s) to have the same number of zeros and poles, which 
explains why the non trivial zeros of -F/v(/Lit; s) do not tend toward — oo as t — ► 1. The fact that 
these zeros tend toward s = seems to support the hypothesis that the evenness of Fjq{p; s) 
stems from the invariance of x + 1/x under the map x i— » 1/x. 

We remark that explicit formulae for fiN(pt] £) and /at (/it; £) may be recovered from H^d^t] s) 
and Fjsr((J,t] s ) via Mellin inversion. In lieu of explicit formulae we give the following qualitative 
corollary Theorem 11.91 This corollary follows immediately from Mellin inversion and we will 
not prove it here. 

Corollary 1.10. For each t e [0,1], fN{nuQ and /i/v (/■**;£) are Laurent polynomials. 
Moreover /jv(jU t ;0 and ir~ N h N (fi t ; arc in Q[£, 

1.5. The Number of Reciprocal Polynomials in X[x] with Bounded Degree and Mahler Measure 

We now turn to an application of the theory of multiplicative distance functions to Dio- 
phantine geometry. The mechanism by which we may infer information about the range of p 
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on 1[x\ from the range of p on R[x] stems from the fact that as T — > oo the cardinality of 
Jj N+1 H TUn is approximately the volume of TUn- Of course we may apply this principle more 
generally, but we limit ourselves to the case of the reciprocal Mahler's measure since reciprocal 
polynomials hold a distinguished role in the study of integer polynomials with small Mahler's 
measure |10j . 

A polynomial of degree M is called reciprocal if f(x) = x f (1 / x) . Each reciprocal polyno- 
mial in Z[x] corresponds to a reciprocal Laurent polynomial in 7L\x + l/x]. We denote the set 
of reciprocal polynomials in 7L\x\ with degree at most N and Mahler's measure less than or 
equal to T by M N {T). 

Theorem 1.11. Let N be a positive integer. Then, as T — > oo, the cardinality of Mn{T) 
satisfies the following asymptotic estimates. 



where the constant implicit in the O-notation is dependent on J. 

Proof. Suppose that / is a reciprocal polynomial in Z[x]. The subset of reciprocal polyno- 
mials ofZ[x] is closed under multiplication. If deg(/) is odd then /(— 1) = 0, and f(x)/(x+l) is 
a reciprocal polynomial of even degree. Furthermore, the multiplicativity of Mahler's measure 
implies 



Thus, when studying the range of values Mahler's measure takes on reciprocal polynomials, it 
suffices to consider only even degrees. We assume that deg(/) — 2 J, and let p{x) — x~ J f(x). 
Clearly p is a reciprocal Laurent polynomial and there exists g(x) G Z[x] such that p{x) = 
g(x + l/x). It follows that = p(g). 

We now turn to Ai n (T) . Notice that M. n (T) consists of polynomials with both even and 
odd degrees. By our previous remarks, if N = 2 J + 1 is odd, then the set of polynomials in 
M n (T) with odd degree is in one to one correspondence with the set 



Likewise the set of polynomials in M n (T) with even degree is in one to one correspondence 



If N — 2 J is even, the set of polynomials in Mn(T) with odd degree is in one to one 
correspondence with the set TUj-i(p) fl Z J_1 , while the set of polynomials in A4n(T) with 
even degree is in one to one correspondence with the set TUj(p) n Z J . 

Then, by well known results, the Lebesgue measure of TUj{p) gives a good approximation 
of the number of integer lattice points contained in TU,j(p) when T is large. Specifically, 





{ 



a e Z J : p(a) < T) = TUj{p) n Z J . 



TUj(p)DZ J . 



# (TUj(p) n Z J ) = \ J+l {U,(p))T J+1 + 0(T J ) as T - oo. 



See [1 Chapter VI, §2] or pQ §12] for details. 



□ 
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This Theorem is useful since we can explicitly compute \j+i(Uj(p)) using Theorem ll.6l For 
example, 



#M (T) =2T + 0(1) 



#M 2 (T) -- 
#Mi{T) = 
#M 6 (T) = 
#M S (T) = 
#Mi (T) 



16^2 
3 



0(T), 
f T 3 + 0(T 2 ), 

ifi§pr 4 + o(T3), 

6^p T 5 + 0(T 4 ); 

_ 2147483648 rpG 
" 21223125 



0(T 5 



#Mi(r)=4T+0(l), 



#-M 3 (T) = 
#-M 5 (T) = 
#-M 7 (T) = 
#-M 9 (T) = 
#M n (T) 



32rp2 _ 

128^3 
5 1 



O(T), 
f 0(T 2 ), 



262144 rp4 
4725 



T 4 + 0(T 3 



1310720 
11907 



T 5 + 0(T 4 ), 



4294967296 ^6 
21223125 



T 6 + ( T 5 



1.6. The Structure of Moment Functions 

The evaluation of Fn(&;s) and Hn($>;s) depends on the multiplicativity of <& as well as 
the specifics of the root function cj>. By exploiting the multiplicativity of $ we may express 
Fn((&;s) and Hn($;s) in fairly simple terms dependent only on 4> and N (and of course s). 
For the remainder of this section we will view <!> as fixed. Many of the structures introduced in 
this section are dependent on <I>, but this dependence will be suppressed in an effort to simplify 
the notation. 

We begin with Hn(s). For each s = a + it with a > N, let rj s be the Borel measure on C 
defined by 

dr, s ( 7 ) = ^( 7 )- 2CT dA 2 (7). 
Next we define a Hermitian form on the Hilbert space L 2 (rj s ) by setting 



(P\Q) 



<t>{ir Zs P{l)Q{l)d\ 2 ( n ) for each P,Q € L z (n s ). 



Notice that when s is real this is just the inner product associated to the norm on L 2 (ij s ). It is 
easy to verify from Theorem 11.11 that the polynomials 1, 7, 7 2 , . . . , j N ~ 1 are in L 2 (rj s ). In fact, 
any complete set of N polynomials, that is a set {Pi, P 2 , . . . , Pn} in C[ 7 ] with deg P n = n — 1, 
is in L 2 (ij s ). 

Theorem 1.12. LetK(s) > N, andlet P = {Pi,P 2 , ■ ■ ■ ,Pn} be any complete set of monic 
polynomials in C[y]. Then, 

H N (^;s) = detWp, 

where Wp is the N x N matrix whose j, k entry is given by Wp[j, k] = (Pj\Pk)- 

The matrix Wp is known as the Gram matrix of the set P with respect to the Hermitian form 
(■|-). When s is real we may view P as spanning a parallelepiped in L 2 (rj s ). As such det Hn(s) 
is the volume of this parallelepiped. Moreover, since Hn(N + 1) is essentially the volume of 
the starbody Vn can also be regarded as the volume of a parallelepiped in the Hilbert space 
L 2 {vn+i)- 

Perhaps the most useful aspect of Theorem 11.121 is that it is independent of the complete 
family of monic polynomials chosen. Thus, a wise choice of P — for instance one which is 
orthogonal with respect to the Hermitian form — may make the evaluation of det Wp easy. Of 
course the coefficients of such orthogonal polynomials will be dependent on s. 
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Corollary 1.13. Let !ft(s) > TV and let Q = {Qi, Q 2 , . . ■ ,Qn} be the complete family of 
monic polynomials specified by 

(Qj\Qk)=^ s (Qk)Skj for j,k=l,...,N. 

Then, 

N 

H N ($;s) = l[m s (Q k ). 

n=l 

When s is real ^t s (Q, i ) is simply the norm squared of Q n in L 2 (rj s ). 

As we shall see the evaluation of Fn(s) is much more complicated, due in part to the fact 
that a polynomial in R[x] may have both real and complex roots. In spite of this difficulty 
results similar to Theorem 11.121 and Corollary 11.131 are available. These can be stated by 
replacing the Hermitian form used in the calculation of Hn(s) with a skew-symmetric bilinear 
form associated to $. The matrix of skew-symmetric bilinear forms formed in analogy with 
Wp is antisymmetric and we will replace the determinant with the Pfaffian — an invariant of 
antisymmetric matrices — in order to give a succinct formulation of Fn(s). 

In analogy with the Hermitian form (-|-) we introduce the skew-symmetric bilinear forms 
(•, -) R and (•, -) c by 



(P,Q)i 

and 



(P,Q) c = -2i 



^x)- s <f>(y)- s P(x)Q(y) sgn(y - x)dxdy, (1.9) 



4>{Pr s 4>(P)- s P(P)Q{P) sgn3(/3) dX 2 (/3), (1.10) 



where as before these bilinear forms are implicitly dependent on s. The skew moniker stems 
from the fact that (Q,P)r = —(P,Q)tr (and similarly for (-,-)c)- When 5R(s) > N it is easily 
verified that the integrals defining (P, Q)r and (P, Q)c converge when P and Q are polynomials 
of degree at most N — 1. We may create another skew-symmetric bilinear form by specifying 
that 

(P,Q) = (P,Q)n + (P,Q)c (1.11) 

Now given any complete family of N monic polynomials P — {Pi, P 2 , . . . , Pn} ^= we may 
create the N x N antisymmetric matrix Up whose j, k entry is given by (Pj,Pk). As before the 
entries of this matrix are functions of s. 

An important invariant of even rank antisymmetric matrices is the Pfaffian. If N — 2J and 
U is an N x N antisymmetric matrix, then the Pfaffian of U is given by 

1 - 

PfC/ =^T7! E sgn(r) 11^(2^-1), t(2j)], (1.12) 

reS N 3 = 1 

where Sn is the symmetric group on {1, 2, ... , N}. The Pfaffian is related to the determinant by 
the formula det U = (Pf U) 2 (see for instance [6l Appendix: Pfaffians]). Thus, psychologically 
at least, the Pfaffian of U may be thought of as the signed square root of the determinant of 
U. 

One of the major results in this manuscript, and the one we will spend the most time 
proving, is that Fn{s) can be represented as the Pfaffian of Up for any complete set P of N 
monic polynomials in Rpy]. However, before this claim can be made it is necessary to adjust 
our definitions for the case when N is odd. 
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Theorem 1.14. Let Sft(s) > N and let J be the integer part of (N + l)/2. If P 
{Pi (7), £2(7), ■ • • i Pn(j)} is any complete set of monic polynomials in C[y] then 

F N ($;s) = PfU P , 

where Up is the 2 J x 2 J antisymmetric matrix whose j, k entry is given by 

{Pj.Pk) ifj,k<N, 



sgn(fc - j) 



(x) s P m in{j,k}( x ) dx otherwise. 



(1.13) 



Notice that when N is even then the first condition in equation (| 1 . 1 3[) always holds. 

As is the case with Hn(s) a smart choice of P yields a simple product formulation for Fn(s). 
Specifically, when N is even we may use a complete family of monic polynomials which are 
skew-orthogonal. 

Corollary 1.15. Suppose that N = 2J, 3?(s) > TV and let Q = {Qx, Q 2 , ■ ■■ , Qn} be any 
complete family of monic polynomials specihed by 

{Q2k-i,Qsj) = -{Qij,Q%k-i) = <h,M.:.(l,i and (Q 2 j,Q 2 k) = {Qij-l,Qik-l) = 0, 
for j,k = 1, . . . , J. Then, 

J 

Fn($;s) = Y[m a (Q j ). 
i=i 

The quantities Wl s (Qj) are referred to as the normalization(s) of Q. 

In the special case of multiplicative distance functions whose root functions satisfy certain 
symmetries we may write -F/v($; s) as a determinant. 

Corollary 1.16. Suppose that U(s) > N and let J be the integer part of (N + l)/2. 
Furthermore suppose that P is a complete family of monic polynomials in M.[x] such that P n 
is even when n — 1 is even, and P n is odd when n — 1 is odd. If the root function of <& satishes 
cj>(-f3) = 4>{f3) for every € C then, 

Fat($;s) = det A P 
where Ap is the J x J matrix whose j, k entry is given by 

A P \j,k] = U P [2j-l,2k}. 



2. The Proof of Theorem [L~I\ 

Since <I> is continuous, non-negative and positive definite, we find that <f> is continuous and 
4>{ct) > for each a G C. The asymptotic properties of root functions are derived from the 
continuity of multiplicative distance functions. To see this, let a and b be nonzero complex 
numbers. By homogeneity, 

<i>(oa; - b) = \a\ $ \ x — - J = \a\ <p 

By continuity lim Q(ax — b) = $(—6) = |6|, and thus 

|a|->0 ' 

lim \a\ <!>(-) =\b\. 
Setting 7 = b/a we see that (^(7) ~ I7I as I7I — > 00. 
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The other direction is more complicated. Suppose that i\> : C — > (0, oo) is a continuous 
function such that ip ~ |-y | as |-y | — ► oo. We will use a modification of Mahler's original proof 
that n is continuous to prove the continuity of the function 

N N 
n— 1 n— 1 

Certainly $ satisfies the other axioms of multiplicative distance functions. 

In fact we will prove that ^ is continuous with respect to the stronger topology induced 
by uniform convergence on compact subsets of C. Suppose that {fk(%)} is a sequence of 
polynomials in C[x] such that 



fk(x) = a kNk Y[ (x - Ikn) for k > 



n=l 

and 



N 



lim f k (x) = f(x) =aY\(x- j n ), 

k — >oc 



n=l 

uniformly on compact subsets of C. We will show that 

lim *(/ fc ) = *(/). 

k — >oo 

By an easy corollary to Hurwitz's Root Theorem (see for instance [8]) we may reorder the 
roots of each fk(x) so that 

lim 7 fcn = 7„ for n = 1, . . . , N. 

k — >oo 

For each k > 0, define the polynomials and h k (x) by, 

g k (x) = a Y[ - Ikn) and h k (x) = (x - j kn ), 

n=l n=N+l 

and notice that g k (x) — > f(x) as k — > oo. Furthermore, since g k (x) is of degree N for all k, it 
follows that this convergence is uniform on compact subsets of C. Now, 

\g k {x)\ \h k (x) - 1| = \f k (x) - g k (x)\ 

< \fk{x)-f(x)\ + \f(x)-g k (x)\, 

from which it follows that {h k (x)} converges to the constant function 1 uniformly on compact 
subsets of C \ {71, . . . , 7jv}. 

The constant coefficient of h k (x) is given by 

a - 1 -■- 

n=N+l 

and thus, by choosing a point x € C\ {71, . . . ,7at} and using the fact that h k (x) — > 1 we have, 

I, 1 1 n=JV+l J 
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This is the key fact needed to prove the theorem. 

lim *(/ fe ) = lim (\a kNk \ TT V(7fen) ) 

k — > oo k — >oo \ I 



\ l. 11 n=JV+l J L n=l 

-^({^J +1 bl " l }{' al S [ 1 * h '" ) } 



Where the last equation is a consequence of the fact that |7fc n | — > oo as fc — > oo, and "0(7) ~ M- 
From jUTJ it follows that, 

{AT ~j N 

m n ^(7^) = m n ^(-rn) = *(/), 
71=1 J 71=1 

where the second equality follows from the continuity of ip and the fact that lim jkn = In ■ 

k — >oo 



3. The Proofs of Theorem IL5I and Theorem IL6I 



We will prove Theorem 11.51 for the real case and leave the complex case to the reader. We 
will view $ as fixed and suppress any notational dependence on $. 



Proof of Theorem 11.51 Let B be the N + 1 dimensional unit ball centered at the origin. 
Then, since Un is bounded we can find a positive constant r\ so that 

U N C r]B and thus TU N C TrjB. 

Let A\/t — {(b, 1/T) : b <= R*}. For instance, .Ai is the hyperplane of coefficient vectors of 
monic polynomials of degree N. It follows that 

(Ax n TU N ) C (Ai n Tr\ B). (3.1) 

Depending on the value of T, the set (A\ PI -B) is either empty or an iV-dimensional ball. It 
follows from (13.11) that 



In{t) < \jf(Ai n Tr]B), 

and there exists an e > such that if T < e then fjy(T) = 0. 

Clearly the set of polynomials with leading coefficient 1/T and distance 1 is given by 
At/T tlU N . Notice that A 1/T = (1/T)A\. Thus we find that 

A 1/T nu N = hA x r\TU N ). 

It is easy to see that (A\/t n Un) — * Un-i as T — > 00. Thus 
\2n{Un-i) = lim Aat(Al/ t n Wat) 

T— *oo 

- lim Aw fi(Anra w )) - lim iv(7 ] 



Again we will prove Theorem II. 61 in the real case and leave the complex case to the reader. 
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Proof of Theorem 11.61 Using the notation from the proof of Theorem 11.51 the volume of 
Un is given by 

X N+1 (U N ) = X N (A T r)U N )dT. 

. R 

By the absolute homogeneity of $ we see 

x N (A T nu N ) = \T\ N X N (Ax n {T^Un) = iri^/jvdTr 1 ). 

And thus 



Xn+i(Mn) 



it v / JV (|r|- 1 )dr = 2 



THfaiT-^dT. 



Finally, by setting £ = T 1 we find 
Xn+i{Mn) = 2 



4. The Proof of Theorem \F7\ 

As remarked previously Theorem 11.71 was originally proved by S-J. Chern and J. Vaaler in 
PQ. Their technique for evaluating Fff(fj,;s) involved a number of rational function identities 
which were specialized to /i. In this section we will present a different proof which relies on 
Theorem II .141 (or rather Corollarv ll.l6j) . 

It is worth remarking that Theorem 11.141 (and its corollaries) reduce the determination 
of Hn($>;s) and Fn(&',s) to the (not necessarily trivial) calculation a number of Hermitian 
forms and skew-symmetric bilinear forms dependent on and then the computation of a 
determinant and a Pfaffian. When $ = /i it is convenient to use the family of monic polynomials 
P = {1, 7, j 2 , . . . , 7 A ' -1 }. For H]y(p; s) the integrals defining the Hermitian forms of pairs of 
elements of P are elementary, and moreover P is orthogonal with respect to the Hermitian 
form. Thus, Hn({i;s) is the determinant of a diagonal matrix with entries that are easily 
computed. The details of this computation are left to the reader (or can be found in [1]). 

The integrals defining the skew-symmetric bilinear forms used in the computation of Fn (/•*; s) 
are slightly more complicated, but still elementary. And since /x satisfies the conditions of 
Corollary 11.161 Fn(/i;s) is given by the determinant of a matrix whose entries are given by 
these skew-symmetric bilinear forms. This matrix is more complicated than the matrix which 
appears in the formulation of ff/v(/x; s) but nonetheless its determinant can be computed. 

Lemma 4.1. Let $ = fx and set J be the integer part of (N + l)/2. Dehne A P to be the 
J x J matrix defined as in Corollarv \1.16l Then, 



1 



A P \j,k] = < 



2k - 2j + 1 \2j-lJ \s-2k 
2 



[ 2j-l \s-2j + l 



if k < 



if k 



N + l 
N + l 



(4.1) 



We will defer the proof of this lemma to see how we may use it to derive the formulation 
for -Fat (/it; s) given in Theorem 1 1.71 Since the second condition defining Ap[j, k] is only realized 
when N is odd, it is sensible to divide the determination of Fjy (/•*; s) into cases depending on 
whether N is even or odd. 
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4.1. The Even N Case 

When N — 2 J, Ap is defined only by the first condition in (|4.ip . We have written this 
suggestively to indicate terms which depend only on the rows or columns of Ap . It follows that 

F N (n; s) = det Ap[ 3 , k] = det B ■ 2 N j JJ — ^ {^fZJ 

where B is the J x J matrix given by B[j, k] = l/(2k — 2j + 1). The matrix B is a Cauchy 
matrix, and using the well-known formula for the determinant of a Cauchy matrix, 

dets=(-i)Go J n (2fc-2jfi /Jnn( 2fc -^ +i )[' ( 4 - 2 ) 

[l<j</c<,7 J/ [j=lfc=l J 

the denominator of which is 

X\\{{2k-2j + l) = {-lp)\ Yl ( 2 (fc_j) + l)(2(fc-i)-l) 

j=l fe=l [l<j</c<J 

=(-i)(=)| n (2(fc-j)+i) 2 ||n(2(j-.?)+i) 

Substituting this into (|4.2[) we find 



And thus. 



which after reindexing yields the formula for F^{^', s) given in Theorem 1 1.71 



4.2. The Odd N Case 

When N is odd we have J = (iV + l)/2. Looking at (|4.ip we may factor out terms dependent 
only on the rows or columns of Ap to write 

F N (n;s) = det A P = 2 N s J jjj {ff^} ' deti ^ ( 4 - 4 ) 

where B' is the J x J matrix given by 



B'\j,k] = 



The determinant of £?' is clearly a rational function of s which we will denote by b(s). Moreover 
as s — > oo we must have b(s) — > 0, from which it follows that b(s) has fewer zeros than poles. 
It is clear from the definition of B' that b(s) has J simple poles located at the positive odd 
integers not exceeding N. It is also easy to see that b(s) has J — 1 zeros located at the positive 



1 



2k -2j + l 
1 

S-27 + 1 



if k < J, 



if k = J. 
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even integers not exceeding N. Thus there exists a constant k such that 

& (-) = - { n } { n - ^) } 

Notice that the zeros of b(s) exactly cancel the poles at even integers which appear in (j4.4|) . 
That is, 

MM = K . 2 »Jn_^(^)J. ( 4.5) 

In order to determine the value of k we must determine b(s) at another value of s, the obvious 
choice being s = 2 J. In this situation 6(2 J) is simply the determinant of the Jx J matrix whose 
j, k entry is given by l/(2k — 2j + 1). That is, 6(2 J) = det-B, the same Cauchy determinant 
that appeared in the even N case. Thus we have 



and by (|P|) 




which after reindexing yields the formulation for -F/v(m; s) given in Theorem 1 1.71 



4.3. The Proof of LemmaWl 



We compute the entries of the matrix Up[j, k] = (j^ 1 , j k x ) under the conditions that j is 
odd and k is even. The root function of /i is 4>{l) — max{l, |7|}, and hence 



max{l, \x\} s max{l,|j/|} s x^ 1 y k 1 sgn(y — x) dx dy 



= 2 



— OO 

y 



<{i,M>- 



^lyiy-V-'y^dxdy. 



But this integral is elementary, since we may divide the domain of integration into regions 
according to where max{l, |x|} and max{l,|y|} are identically one. The integrals converge 
when SR(s) > j + k. Putting the result into partial fractions form (as a function of s) we find, 



,1 



2s -j-k \ j 



]{j + k) 



2k 



k \j{k-j) 



(4.6) 
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Now, 



= -2i 
= -U 



max 



{lA(3\}- 2s (Py- 1 P k - 1 sgn^(f3)d\ 2 (P) 



max{l,r} 2 V+ fc 1 dr x 
max{l,r}~ 2 V +fe ~ 1 dr x 



*7T 


r 27T-l 






. 


7T J 



The integrals in this expression are elementary, and when 3?(s) > j + fc, we find 

(7^ 1 ,7 fc " 1 )c = -: 



2s-j-k\j-kJ (k-j)(k + j) 
Notice that the first term in (7 J ' — 1 ,7 )c exactly cancels O in (|4.6p . That is, 
= /«i-i«*-l\= 4 2 / 2fc 



^p[j^] = ( 7 -'- 1 , 7 fc - 1 ) = — 



i(j + fc) s - fc \j(k - j) ) (k- j)(k + j) 
4 (a 



j(j -k) \k- s / 

When k < N, the entries of Ap[j, k] are given by Up[2j — 1, 2k]. When k = N + 1 (which can 
only occur when J is odd) we have 

2s 



A P \j,N+l] = 



max{l, |x|} s x 2] 2 dx — 



(2j-l)(s-2j + l)' 



which finishes the proof of the lemma. 

5 . The Proof of Theorem PI 



Our strategy is the same as in the proof of Theorem 11.71 First compute the Hermitian and 
skew-symmetric bilinear forms for a complete family of polynomials and then compute the 
determinant and Pfaffian of the appropriate matrices whose entries are these bilinear forms. 
To compute Fm{p\ s) we will use Corollarv l 1 . 1 61 with P = {1, 7, 7 2 , . . . , j N ~ 1 }. The formulation 
of Hm{p\ s) given in Theorem ll.Sl will not be presented here since it (or rather a minor variation 
of it) is the subject of [9]. 

Lemma 5.1. Let $ = p and let J be the integer part of (N + l)/2. Define A P to be the 
J x J matrix defined as in Corollarv l 1.1b] Then, if k < (N + l)/2, 



2k -I 

k — n 



A P [j,k] = Y, 

71=1 

and, 

A P [j,^] = f 
where 



16s 2 



s 2 - (2n) 2 



E 

m— 1 



n=l 



N 

JV+l 



2 N-2 



E 



2j-2' 

j — m 

'2j-2 

j - m 



2n 



2m -I \{2n) 2 - (2m- l) 2 



2n 



1 



2m- 1 \s 2 - {2m -If 



M 
rn 



M 

TO — 1 



, (5.1) 



(5.2) 



(5.3) 



Lemma 15.11 is proved by brute force, and in fact much of the proof involves massaging the 
entries of Ap into the form given in the statement of the lemma. This form is not the most 
natural but will be useful for our purposes as we ultimately need to take the determinant of 
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Ap . We defer the proof of Lemma 15.11 in order to see how it may be used to compute the 
formulation of -F/v(p; s) given in the statement of Theorem II. 81 
We define the J x J matrices B, C and D by 



B[m, n] = 



2n 



2m - 1 
2n 



1 



(2n) 2 - (2m - If 



if n < 



N + l 
2 ' 



1 



2j-2 
j -m 



2m - 1 \s 2 - (2m- l) 2 



and Ji] = 



if n = ^l, 





'2k- 


- 1" 


< 


k - 


n 




'2k- 


- 1" 


< 


k- 


n 



16s 2 



s 2 - {2nf 



s 

4> 



if fc<^± 



if fc = ^±l. 



Thus Ap = DCB and Fn(p; s) = det Ap = det -B • det C • det D. This is convenient since, the 



2j-2 
3— m 



= 0). 



matrices C and £> are triangular (since for instance if m > j then 

When iV is even we have J < (N + l)/2 and hence the first conditions defining B and D 
hold. Thus det B is simply a rational number, and thus computing the diagonal entries of C 
and D we see that there is some rational number vn so that 



Fn(p\s) 



j 

VN Ii 

3 = 1 



(2j) 2 



J 

vnH 

j 







{N-2jf 



(5.4) 



When N is odd, then the determinant of B is a rational function of s, which we will denote 
b(s). From the definition of B it is easily seen that b(s) is an even rational function with simple 
poles at the integers ±1, ±3, . . . , ±JV. Moreover when s = ±2, ±4, . . . , ±(N — 1) the matrix 
B is singular and hence b(s) = for these values of s. Also from the definition of B it is 
seen that b(s) — > as s — * oo. We conclude that there b(s) has fewer zeros than poles. We 
have identified all the poles of b(s), and since b(s) is even we have also identified the complete 
list of zeros of b(s). Notice that the zeros of b(s) exactly cancel the poles which arise from 
the diagonal entries of D. Putting these observations together we find that there exists some 
rational number (which we also denote vn) such that 



Fn(p;s) 



j 

VN Ii 

3= 



(2j-l) : 



= v N 



n 



3- 



(N - 2jf 



(5.5) 



We may find the value of by explicitly computing the determinant of B, C and D, and 
noting that det-B is a Cauchy determinant. However, by casting the constant in another 
context we may find its value in the literature. 

Equation (|1.8[) implies that vn — \n(Un(p))- That is vn is the volume of the set of b G M. N 
with jo(b) = 1. This observation is useful since the volume of Um{p) has been computed by 
S. DiPippo and E. Howe in pH Proposition 2.2.1]. DiPippo and Howe show that 



Xn{U n {p)) 



N N 



2n 



2n-l 



JV+l- 



(5.6) 



In fact, DiPippo and Howe report this number as the volume of the set of monic coefficient 
vectors of polynomials of degree N with real coefficients and all roots on the unit circle, however 
in the^ course of their computation they show that this volume exactly equals the volume of 
Ajv(Wjv(p)). Putting (|5.6[) together with (|5.4p and (|5.5p we arrive at the formulation of Fjy(p; s) 
given in the statement of Theorem 11.81 
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5.1. The Proof of Lemma I5.il 

Throughout this section we set $ = p. 

Before proving Lemma 15.11 a few results about the binomial-like coefficients are in order. 
First let us see how these coefficients come about. The entries of Ap are defined by integrals 
in which factors like <^(7)~ S 7 J occur. In order to evaluate these integrals it is convenient to 
use the change of variables 7 1— > 7 + I/7 since ^(7 + I/7) = max{|7|, ^l" 1 }. We are left with 
integrands of the form 



7 



3-1 



1 



r 



This is beneficial since we may use the binomial theorem to expand the latter as a finite sum. 
It is this expansion together with the Jacobian of the change of variables which produce the 
variants of binomial coefficients given in (|5.3|) into our calculations. A few facts regarding these 
coefficients are necessary. 

Lemma 5.2. Let j and k be positive integers. Then, 



(i) [x + 

(ii) 2* = £ 



71=0 





1 








X 


'k-1 


(k 


n 



E 



m 



x 



j-2m 



2 1 \ — • 

(iii) If j is odd, then — = y 



J 



771 = 



(iv) If j is odd and k is even, then 
-k i fc 



2 J ' 



j(j + k) 



EE 



1 




1 




I 


j " 


- 2m' 


3 


- 1" 




"fe- 1" 


m 




n 



1 



1 



j — 2m j — 2m + k — 2n 



Proof. To prove (i) we use the Binomial Theorem to expand (x + 1/ x) 3 1 (x — 1/x) and 
collect together terms with like powers of x. Fact (ii) follows by taking the derivative of both 



sides of (i) and setting x = 1. 



consider the path integral 



To prove (iii) let w be a path in the complex plane that does not pass through z = and 

(5.7) 



1 



l-—)dz. 



If j is odd then (i) implies that the integrand consists of even powers of x. Thus the integral 
in 15.71 depends only on the end points of w. Notice then that, 



9.7 
j 



x-i 1 dx = — 
2 







I - 


-I) 


Jul \ 





3-1 



1_ ^2 ) "' ; 



where the second equality follows from the change of variables x 1— > z + 1/z, and lu is any path 
in the complex plane starting at z = — 1 ending at z = 1 and not passing through z = 0. Using 
(i) and the Fundamental Theorem of Calculus we find 
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To prove (iv) we notice that 

23 +k 



j(J + k) 



.k-l 



c- 7 1 dx dy 



1 i 



"i-i" 




r 2 


■<t>+(y) 










m 




-2 


-i 



(J- 2 ™- 1 dxdy, 



where the second equality stems from the change of variables x i— > x + 1/x, and 4>+(y) = 
(y + vT-4)/2. Again we use the fact that j is odd to conclude that the resulting integral is 
path independent. Assuming that k is even we may evaluate the inner integral and simplify to 
find 



23 +k 

j(j + k) 



- V 

2 ^ 

m—0 
3 k 

rn—0 n—0 
m—0 n—0 



1" 


1 


r 2 




I 


j ~ 




2m . 


-2 


J 


- 1" 




'fo- 


1" 


l 


m 




ri 




j — 2m 


'j 


- 1" 




'k- 


1" 


1 


m 




n 




j — 2m 



_/'+fc-2m-2n-l 



d;y 



Proof of Lemma POl We will use brute force to compute the entries of the matrix Up [j, k] = 
(7 J_1 ,7 fe_1 ) under the conditions that j is odd and k is even. 
We will first evaluate (7 J ' _1 ,7 fc_1 )R. Define the functions </>-, 0+ : K — * C by 



0_(a) 



a — Va 2 — 4 



and + (a) 



where is a fixed branch of the square root which maps the positive real axis to itself. The 
root function of p restricted to the real axis can be given by, 

-0_(a) if a < -2, 
(a) = { 1 if-2<a<2, 
</>+(a) if a > 2. 



From the definition of (7-? 1 ,7 fc 1 )r it follows that 



</>(y) s y 



c) s a;-* 1 da; dy. 



Hv) 



The inner integral of (|5.8p can be written as 



(-^_(x))-*aB* _1 di 



c) x J 1 dx = < 



^(-2) 



^(2) 



0+ (x) s x J 1 dx 



y< -2, 
-2 < y < 2, 

2 < y. 



(5.8) 
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Each of these integrals converge when $R(s) > j. When y < 2 we may use the change of variables 
x i— > x + 1/x to write 



Hv) 



<t>-(v) 



— OO 
J 



(x-\ 






dx 




--)'" 


H) 






X 




X 



E 

m=0 



1 r^-(y) 



where the second equation comes from Lemma 15.21 Similarly, when y > 2 we may write 



(s)-V- 1 da; = 51 



m=0 



i-i 

m 



- Sx j-2m-l dx _ 



Evaluating these integrals using the fact that j is odd we find, 

' ' {-^-(l/)}-V-(») j - am 



E 

m=0 



3-1 

m 



Hv) = { 



y J + 2i 
j 



E 



m=0 



— s + j — 2m 



m 



V. J m =0 



i-i 

rn 



—s + j — 2m 
— s + j — 2m 



2/<-2, 



-2 < y < 2, (5.9) 



2<y. 



Now (7 J 1 ,7 fe 1 )r is given by 

{-4>-(y)r s y k - 1 F(y)dy + 



y k - i T(y)dy + 



ct> + {y)- s y k - l Hy)dy). (5.10) 



Using (|5.9p we find 

» = E 



i-i 

rn 



1 



{-0-(»)}- 3 "0-(v) J '- 2m W fc - 1 d|/. 



— s + j — 2m 

The change of variables y i— > y + 1 /y together with Lemma 15.21 yields 
3 k r 



« = EE 



'j - 1" 




"fc- 1" 


m 




n 



1 



— s + j ; — 2m 



{-y}-^V + dy. (5.11) 



m— n— 

Again substituting (|5.9|) into (|5. 10|) we may write the © as 



2 .7 



E 

m=0 



i-i 

m 



1 



— s +j— 2m 



V k 1 dy + - 

-2 J 



since k is even 



y* +fc - 1 d» (5.12) 
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Similarly © can be written as 



m— C 
3 

+ EE 



.7-1 
m 



1 k 

- — y 

n — 2m ^ — ' 



-s + j — 2m 



k - 1 

n 



y 



-s+k-2n-l 



dy 



m— n— 

2^' +1 

J n=0 



"i-i" 




"fc- 1" 


m 




71 



1 



— s + j — 2m 
s+fe - 2 "- 1 dy. 



-2s+j-2m+k-2n-l 



dy 



Evaluating the integrals in (|5.1ip . (|5.12|) and (|5. 13[) we may write (-f j 1 , / y k 1 )r as 



2 J+2 k 



~3 - 1" 




~k-l~ 


m 




n 



2J + fe+2 



— E 

n=0 



k - 1 

n 



1 



j k 



s — k + 2n 



4 EE 

m— n— 



- j - 


k 


+ 2m - 


"i-i" 




~k - 1" 


m 




n 



1 



1 
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(5.13) 



s — j + 2m s — fc — 2n 



Decomposing © + © into partial fractions gives 

k 

E 



+4EE 

//7=0 n— 
j A: 

-sEE 



2 J - 


-fe 


+2 


iC? ■+ 


*) + " 


"i-i" 




"fc - 1" 


m 




n 






"fc - 1" 


m 




n 



fe-1 

n 



1 



1 



s — A: + 2n 
1 



2m — 2n + k — j s — fc + 2n 
1 1 



2m -2n + k- j 2s - k - j + 2n + 2m 



© 



We now turn our attention to (7^ ,7 ) 



-2/ 



^ 1 (/3)- 2s (^- 1 /3 fc - 1 sgn3(/3) dA ( /3). 



After the change of variables /3 1— ► /3 + 1//3 we may rewrite (7 3 1 , 7 \ as 



{|/3|,|/3- 1 |}- 2s (/3 + J 



fc-i 



x sgn Sy I /3 + 



1\ d\ 2 (J3) 







(5.14) 



The integrand is invariant under the map (3 1— ► 1//?, and thus we may replace the domain of 
integration with C\D (recall that D is the closed unit disk). In this domain, sgn(3(/3+l/ /?)) = 1 
if (3 in the open upper half plane, and is equal to —1 if (3 is in the open lower half plane. After 
an easy simplification we may use these facts to rewrite (7 J_ ,7 )c as 



-4/ 



H\D 



3-1 



n*~V, 1 



' 1\ dA 2 Q3) 
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where H is the open upper half plane. Employing Lemma 15.21 we may rewrite this as 

j k 

m— On- H\D 



"i-i" 








m 









Switching to polar coordinates this becomes 

~ 4j ££ J 1 * 1 e( 2m - 2n+fc -^ ifl dfl 

rn — n— 



"i-i" 




"fc- 1" 




m 




n 





-2s+k+j-2n-2m-l ^ r 



Of course, these integrals are elementary and we finally can write (7^ 1 ,"f k 1 )c as 

1 1 



EE 



"i-i" 




"fc - 1" 


m 




n 



m—0 n—0 

Notice that this exactly cancels © from (|5.14[) . Thus, 

© 



2m - 2n + k - j 2s - k - j + 2n + 2m 




j(j + k ) J 



n=0 



1 



s — k + 2n 



j k 

+4££ 

rn— n—0 



"i-i" 




k-1 


m 




n 



1 



1 



2m — 2n + k — j s — k + 2n 
Using Lemma IS~^1 we may replace and © so that 

1 1 



<7 J - 1 ,7 fc - 1 >=4^^ 

m—0 n—0 

1 



j'-i 

rn 



k - 1 
n 



1 



j — 2m j — 2m + k — 2n 
1 1 



x 



j — 2m s — k + 27i 2to — 2n + k — j s — k + 2n 

Reindexing by m h-> j — 2m and n 1— > fc — 2n allows us to write (7 J ' — 1 , 7 fc l ) as 
i fc 



i-i 



m=—j n— —k L 2 



Notice that since j is odd, 



fc-1 

k — n 
2 



3-1 
372 



1 



1 



m(m + n) m(s — n) (n — m)(s — n) 



(5.15) 



0, so we need not worry about the denominator of the 



summand being identically zero. Next we use the fact that, for instance, 



J-l" 




"i-i" 


j+m 




3— m 


L 2 - 




1- 2 J 



to index the sums in (|5.15|) over only positive integers. Doing this and simplifying the resulting 
summand we find, 



k j 



i-i 

j-m 
2 



fc-1 

k — n 
2 



1G 



The summand is identically zero if m > j or n > A; thus we may replace the upper bounds 
of summation with 2J. Likewise the summand is identically zero unless n is even and m is 
odd. We may thus may reindex the sums by m <—> 2m — 1 and n 1— > 2n. Making these changes 
and simplifying the resulting expression we arrive at the formulation of Ap\j, k] given in the 
statement of the lemma in the case where fc < (N + l)/2. 
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When k = (N + l)/2 we have 



25 



"'-J -1 da 



(a) a 



"at 3 ' 1 da + 



r 7 1 <ia 



+ (a) s oP 1 da 



2 j+i 
© 



"2£ 



i-i 

m 



- s+j -2m-l da ^ 



by the same change of variables used before. Replacing © with the formula given in Lemma 
15.21 and reindexing the sum by m i— > j — 2m we find 



m=-j 
3 



i-i 

2 

i-i 

j—m 
2 



1 



1 



s — m to 



s — m s + to to 



Since j is assumed to be odd we may make the substitution j i— > 2j — 1. Similarly, since the 
summand is only non-zero when to is odd we may reindex the sum bymH 2to — 1. Thus, 
after simplifying, 



2j-2 

j — TO 



1 



2m - 1 \s 2 - (2m - 1) 



where the change in the upper index of summation is justified since 



Now from (ii) of Lemma l5.2l we have 

2J r 



n=0 



2J r 



and thus, since J = (N + l)/2 we have 



N+l 



n=0 



iV 

N+l _ 



n=0 



2n. 



2J- 1 
J — n 



2j-2 
j-m 



2n, 



if m > j. 



It follows that 

Ap[j,W=2- w + 2 



N 

N+l _ 



m— 1 



2j-2 
j-m 



1 



2m- 1 \s 2 - (2m - 1) 



which after reorganization yields the formula for Ap[j, (N + l)/2] given in the statement of the 
lemma. □ 



6 . The Proofs of Theorems Ol and 01 



The proofs of Theorem 11.121 and especially Theorem 11.141 are rather technical. In order to 
see past the technical details it is worthwhile to look at the general strategy for these proofs. 
Looking at Fn($>; s) and Hjy(Q; s) we see both integrals are of the form 

$(b)" 2s dA(b), 



monic 
coefficients 
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where, of course, the monic coefficient vectors we are integrating over and the measure A 
are dependent on whether we are looking at the real or complex moment function. In order 
to evaluate this integral we need to exploit the multiplicativity of $ by making a change of 
variables which allows us to integrate over the roots of monic polynomials as opposed to the 
coefficients. That is, we use maps of the sort 

E : roots —* monic coefficients, 

to write something of the form 



N 



n 



JacE(7) dA'(7), 



roots 



= 1 



where A' is the appropriate measure on the space of roots. At this point we begin to see 
difficulties arising in the case of real moment functions which do not occur for complex moment 
functions. Namely, the space of roots of real polynomials of degree N is more complicated then 
the space of roots of complex polynomials of degree N. To be quite explicit, the space of roots 
of complex polynomials of degree iV is essentially just the identification space formed from the 
canonical action of Sm on C^. Consequently, 

' N ~| 

n<K7n)- 2s Mac.E(7) d\ 2N ( 7 ). 

vn=l J 

The space of roots of real polynomials of degree N on the other hand is partitioned into 
components determined by the possible numbers of real and complex conjugate pairs of roots. 
That is, 



C N 



L t M>0 
L + 2M = N 



where -Fl.m^; s) is given by 
1 



2 M L\M\ 




M 



Jl <M/3)- s <M/?r s \ Ja,cE(<x,p)d\ L ((x)d\ 2N ((3), 



m—l 



That is, Fl m (^i s ) measures the contribution to Fm(&; s) of polynomials with L real roots and 
M complex conjugate pairs of roots. In the case of Hjy(&; s) we have chosen to integrate over 
all root vectors in C N instead of C N /S N . The 1/N\ term in front of the integral compensates 
for the fact that almost every polynomial gets counted AH times by doing this. Similarly the 
1/(2 M L!M!) allows us to integrate over vectors of roots in the expression for Fm{&; s). 

At this point it makes sense to resolve the ambiguity with the maps represented by E: Let 
En to denote the change of variables from vectors of complex roots to vectors of complex 
monic coefficient vectors, and let El,m to denote the change of variables from vectors of L real 
roots and M pairs of complex conjugate pairs of roots to real monic coefficient vectors. The 
Jacobians of Em and -El.m are related to the Vandermonde determinant. To be explicit, given 
7 G C N let V* to be the N x N matrix whose j, k entry is given by k] — 7j 1 (that is 

V" 1 is the N x N Vandermonde matrix in the variables 71,72, ■ • ■ ,7jv)- It is well known that 
the Jacobian of Em at 7 is given by |dety ) '| 2 . Perhaps less well known, and the content of 
Lemma 17^1 is that the Jacobian of £z,,m a t (ot, (3) is given by 2 M | det y a '^| where we interpret 
(a, (3) as the vector (/3i, . . . , Pm, Pm, «i, ■ ■ ■ , "l, )■ Using this we may write Hm{&] s) and 
F l ,m($;s) as 

H N ($;s) = |l]>fr™r 2S } Idetn 2 d\ 2N (l), (6-1) 
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and 

L M ~\ 

II^r S II <t>{P)~*<t>(PT a \detV a ^\ d\ L ( a )dX 2M ((3). 

£=1 m=l ) 

(6.2) 

Here we can see another source of complexity in the evaluation of Fn($>; s) which does not arise 
in the evaluation of Hn(Q; s). Namely, in the expression for Fl,m(&', s) we have the absolute 
value of a Vandermonde determinant, while in the expression for Hn(<&; s) we have the modulus 
squared of a Vandermonde determinant. In the latter case, we may treat | det V~t\ 2 uniformly 
for each 7 € C N by writing 

I det y | 2 = dctV -detVT (6.3) 

We do not have this luxury when working with Fl ,m($; s). For every (a,f3) E K L x C M , 
det V a '@ is a complex number which is either real or purely imaginary, and we must treat 
each (a, (3) differently based on whether detV a 'P is real and positive, real and negative, on 
the positive imaginary axis or on the negative imaginary axis. It is at this key point that our 
evaluation of Fn(&; s) diverges from the evaluation of -Fjv(/x; s) given by Chern and Vaaler. 



1 



L\M\ 



6.1. The Proof of Theorem [77121 

The evaluation of Hjy(&; s) will give us insight into the evaluation of Fl,m(®', s). So far we 
have only used the multiplicativity of $ to produce a product over the roots of a polynomial. 
To exploit the appearance of this product we will expand both determinants on the right hand 
side of (|6.3[) as a sum over Sjy, 

|detn 2 = f E ^n(a)f[^-A( £ sgn(r) f[ ^-1 \ 

\ o-eSjv n=l ) \reS N n=l ' 



dot V~t dot V 



Substituting this expression into (|6.ip . using the linearity of the integral and combining the 
products gives 



1 r f N 

Hn($;s) = — J] J] sgn((r) sgn(r) ^ I ^ 



Using Fubini's Theorem we may finally see the full usefulness of the multiplicativity of $, 

1 N 
H N (<P;s) = — E sgn(a)sgn(r) n<7 CT( " ) " 1 |7 T(n) - 1 ). (6.4) 

Of course, we must justify the use of Fubini's Theorem, but as remarked previously since 
yr(7i)-i an( j yr(n)-i are polynomials f degree less that iV, they are in L 2 {v s ) for 5ft(s) > iV. 
Now, by reindexing the product by n 1— > r _1 (n) and noting that sgn(u) sgn(r) = sgn(<r o r _1 ) 
we see the right hand side of (|6.4p is simply the determinant of the N x N matrix whose 
j, k entry is given by (7 J_1 ,7 ) (see [9l Lemma 3.1] for details). That is, we have proved 
Theorem 11.121 in the special case where P = {1, 7, 7 2 , . . . , 7 Ar ~ 1 }- In fact, the general case 
is trivially different from this case by noticing that if P = {Pi (7), ^2(7), ■ ■ ■ , P/v(7)} is any 
complete set of monic polynomials and V P '' y is the N x N matrix whose j, k entry is given by 
V p ^[j, k] = P fe (7j) then det V F " = det V . 
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6.2. Remarks on the Proof of Theorem ] 1.1 4\ 

Returning to Fl,,m{®'i s) let us see how the Pfaffian arises in the formulation of Fjv($; s). To 
do this we will use the familiar formula for the Vandermonde determinant given by 

tetV> = 11 (7»-7m)- (6.5) 

l<m<n<N 

Setting 7 = (a, (3) — (/3i, Pi, . . . , (3m, (3m, ol\,..., ccl) we can use (|6.5[) to determine whether 
detV a 'P is on the positive real axis, negative real axis, positive imaginary axis or negative 
imaginary axis. In Lemma l7.3l we will demonstrate that 

|detF a ^| = (-i) M J J] sgn(a fc - aj ) Y[ sgn3(/3 m ) 1 detV a ' . 

\l<j<k<L m=l J 

Next we introduce the L x L antisymmetric matrix T a whose j, k entry is given by T a \j, k] = 
sgn(afc — aj). When L is even, we have the very important identity ( Lemma 17. 4[) 

II sgnK-oy) = PfT a . 

l<j<k<L 

When N is even so is L, and thus in this situation, 

|detF a ^| = {-i) M j \{ sgn3(/3 m )|pfT a -det^ a ^. (6.6) 

In the case when TV (and hence L) is odd we will need to modify our approach since the Pfaffian 
is only defined for even by even square antisymmetric matrices. The point of this section is to 
see the general mechanism which produces the Pfaffian structure in i 7 V(4>; s); for now we will 
assume that we are in the easier case where N is even. 
Substituting (|6.6p into (|6.2p allows us to write 

(-i\ M r f L M 1 

flm®-, s) = n n m- s m- s 

gt xC M U=l m=l ) 

x Pf T a ■ det y dA L (a) d\ 2M {(3) ■ 

At first glance this does not look to be much of an improvement over (|6.2p . however the Pfaffian 
admits an expansion as a sum over Sn similar to that of the determinant (equation 1 1.12p . This 
expansion together with the Laplace expansion of det V a '@ (using minors which depend only 
on a or (3) will allow us to separate this expression into a sum over complementary minors 
whose summand is the product of two integrals: one over K L and the other over C M . The 
integrals over M. L will evaluate to Pfaffians of L x L antisymmetric matrices, the entries of 
which are skew-symmetric bilinear forms of the form (a 3- , a )r. Similarly the integrals 
over C M evaluate to Pfaffians of M x M antisymmetric matrices with entries of the form 
(/3 J_1 ,/3 )c< The evaluation of these integrals is again dependent on Fubini's Theorem. The 
combinatorics necessary to reduce the resulting sum of products of Pfaffians of antisymmetric 
matrices to the formula given in Theorem 11.141 is achieved using a combinatorial formula for 
the Pfaffian of a sum of antisymmetric matrices. That is, the Pfaffian of a sum can be written 
as the sum of a product of Pfaffians. 



7. The Proof of Theorem U~M 

The proof of Theorem [TTTJ] relies on several technical lemmas. In order to clearly see the chain 
of reasoning used to prove Theorem 11.141 we will defer the proofs of these technical lemmas 
until later. First though, we must introduce some definitions and notation. 



Lt(n) 
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7.1. Definitions and Notation 

For each K < N we define 3^ to be the set of increasing functions from {1,2,..., K} to 
{1,2..., A^}. That is, 

3£ = {{l,2,...,A-}- L >{l,2,...,JV} : t(l) <t(2) < ••• <t(JQ}- 

Associated to each te 3^ there exists a unique t' G 3jv-if sucn ^ na * the i ma g es of t and t' are 
disjoint. Each t G 3^ induces a unique permutation t t G Sjv given by 

i(n) if 1 < n < if, 

t'(n-A") if K < n < N. 

We define the sign of t by setting sgn(t) = sgn(tt). The identity map in 3% is denoted by i. To 
each t we associate the subset of the symmetric group given by 

Sjv(t) ={t£ Sn ■ r(k) is in the image of t for k = 1, 2, . . . , AT}. 

For each r € 5jv(t) define the permutations o> G Sk and 7r T G Sn-k by specifying that 

cr T (fc) = r 1 (T(jfe)) fc = i,2,...,A: 

tt t (£) = t'- 1 (r(A' + £)) £ = 1, 2, . . . , A^ - K. 
We may use these definitions to give an alternative (and useful) description of the sign of 

Lemma 7.1. For every t G 3^ and r G Sjv(t), 

sgn(r) 



sgn(t) 



sgn(er r )sgn(7r T )' 



Proof. Clearly i t G S N {t), and i t _1 o t permutes {1, 2, . . . , K} and {K + 1, AT + 2, . . . , A^} 
disjointly. The action of this permutation on {1, 2, . . . , K} is exactly that given by <r T . Similarly, 

ir T (£) = (l^ 1 o T )(K + £) - K for £=1,2,...,M. 

It follows that the cycles in the cycle decomposition of i^ 1 or are in one-to-one correspon- 
dence with the cycles in the cycle decomposition of <r T together with the cycles in the cycle 
decomposition of n T . This yeilds, 

sgnOr 1 o t) = sgn(cr r ) sgn(7r r ). 

In other words, 

sgn(r) 



sgn(t) = sgn(tt) 



sgn(cr r ) sgn(?r T ) ' 



Given an JV x iV matrix W and u, t G 3^, define W u ,t to be the K x K minor whose j, k 
entry is given by W Ui t[j, k] — W[u(j),i(k)]. The complimentary minor is given by W U ',t'- As 
an example of the utility of this notation, the Laplace expansion of the determinant can be 
written as 

det W = sgn(u) ^2 s S n ( t ) det ^u.t ■ det W u , >v , (7.1) 

where u is any fixed element of 3%. We will also use the abbreviated notation W u for W UiU ; 
this is useful notation for working with Pfaffians since if W is an antisymmetric matrix then 
minors of the form W u are also antisymmetric. 

Throughout this section L and M will be non-negative integers such that L + 2M = N. We 
also set P to be a fixed complete family of monic polynomials. We will reserve J for the integer 
part of (N + l)/2, and we will set K to the integer part of (L + l)/2 so that 2AT + 2M = 2 J. 
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We will use a £ R L for a vector of real variables and j3 S C M for a vector of non-real 
complex variables. As before, V a '@ will represent the N x N Vandermonde matrix in the 
variables . . . ,0m, 0M, oh, . . . , Ol and El.m : K L x C M — > M. N is the map given by 

El,m(&,0) — b where 

N L M 

x N + Y, b n x N - n = Y[(x - a t ) n (as - /3 m )(z - ^). 

n— 1 ^—1 m— 1 

It is easily seen that almost every b € R w corresponds to 2 M M ILl preimages under the map 
El,m- 

Matrices will be denoted by capital roman letters, subscripts will be used to define minors 
of a matrix, while superscripts will be used to reflect any variables or parameters on which the 
entries of the matrix are dependent. Thus, for instance W°[ is a minor of W with entries that 
depend on a and (3. 

There are complications in the proof of Theorem 11.141 for odd N which do not arise in the 
even N case. In spite of this disparity we will present the even and odd cases simultaneously. 
Any structures necessary for the odd N case but unnecessary for the even N case will be 
subscripted by o. 

7.2. Steps in the Proof 

As suggested in Section |5] we will use the change(s) of variables El,m- Let 2?l,m represent 
the subset of M. N which consists of coefficient vectors of monic polynomials of degree N with 
L real roots and M pairs of (non-real) complex conjugate roots. That is, T>l,m is the image in 
K.^ of El,m- Clearly is the disjoint union of T>l,m over all pairs of non-negative integers 
with L + 2M = N. Thus, 

F N ($;s) = 

{L,M) JT>L - M 

L \ ( M 




II { IT 0(Pm)- S mn)- S (7.2) 



x J&c E LtM {ct,f3) d\ L (a)d\ 2 M(fl), 

where the sum over (L, M) is understood to be over all non-negative integers L and M such 
that L + 2M = N. 

Lemma 7.2. The Jacobian of El,m js given by 

Jac£ L:M (a,/3) = 2 M jdetV^I . 

Lemma 7.3. Let 7 e C N be given by 

7 = ■ . ■ ,0m,0m,oli, . . . ,ol), 

and let W a 'P be the N x N matrix whose j, k entry is given by 

W a << 3 [j,k} = P k ( lj ). 
Then, if i S js t ne identity map on {1,2, ... , 2M}, 

|detV^| = sgnW j det WftH) M II sg^O 3 ™)} { detWfo J] sgn(a fc - aA 

ted", I m=l ) { j<k ) 
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where as suggested by the notation, the minors Wf. [ and W^ t i , ofW a ^ are dependent only on 
f3 and a respectively. 



Using Lemma 17.21 and Lemma 17.31 we may rewrite (|7.2p as 

*m*;-)= e ^ e f f (n^rUn 

(l,m) ' 'tea", ol /m ) U=i 



M 



x I det W x ? v J] sgn(a fc - a,) U det W? t (-i) M [] sgn3(/3 m ) \ d\ L (a) dX 2M (f3), 

{ j<k J I m=l J 

and Fubini's Theorem yields 

F N @; s) = £ E s S n (t)^ f | II <^)~ s II s g n ^ ~ a ^ \ det rfA i(«) 

| II 0(/3m)^ s 0(^) _s sgn9(/3 m ) I detWg dA 2M (/3). (7.3) 



M! 



v m— 1 



Lemma 7.4. Let K be the integer part of (L + l)/2. Define T a to be the 2K x 2K 
antisymmetric matrix whose j, k entry is given by 



Then, 



Proof. See [2]. 



T a \j fel = <! sgn ( afc ~~ a ^ if 3,k < L + l, 
1 sgn(fc — j) otherwise. 



II sgnK =P1T C 

l<j<fc<L 



□ 



It is worth remarking that when L is even, the first condition defining T a is always in force. 
Since the Pfaffian is only defined for even rank antisymmetric matrices, the second condition 
is used when L is odd to create a 2K x 2K antisymmetric matrix from an L x L matrix. 

Using Lemma 17.41 we may rewrite (|7.3p as 

/ T \ 

1 



Ml 



M 



(l,m) tea~ M 



J] <X/3 m )- s <M/3 m )- s sgn3(/3 m ) } detWg dA 2A f(/3). (7.4) 



It is necessary for our calculations to replace the t <E J^m with elements of J|m • Each t S 3^/ 
induces a unique t £ J^m ^y setting t = t . Notice that t' and t' a differ in the fact that if 
N =/= 2J then t' D (2J - 2M) = 2 J. Clearly sgn(t ) = sgn(t). 

Lemma 7.5. Let R be the 2 J x 2J matrix whose j, k entry is given by 

(P 3 ,P k }m if j,k<N+l 



R\j, k] 



sgn(fc - j) 



(a) s P m in{j.k}(o:) da otherwise, 



32 CHRISTOPHER D. SINCLAIR 

and suppose that t S 3^1 ■ Then, 



1 



L 



J\(t>{u t )- S \ PfT" detW t ? t , d\ L (a) = Pfi? t 



1=1 



When N is odd and t € then Re is an odd by odd matrix. The introduction of t is 
useful since the Pfaffian of Re is defined. 

o 

Lemma 7.6. Let C be the 2 J x 2 J matrix whose j, k entry is given by 

' (Pj,P k )c if j,k<N+l 



C[j, k] = 

and suppose that t S ■ Then, 



otherwise, 



H) 



M 



M\ 

c 



M 



11 0(/3 m )-^(/3 m )- s sgn3(/3 m ) I detwf t d\ 2M (P) = Pf Ct, 



Using Lemma 17.51 and Lemma 17.61 we may rewrite (|7.4[) as 

F N (<f>;s) = E sgn(to)Pf^ PfC to . (7.5) 

(l.m) tel» u 

If u € 3\ J M then either 2 J is in the image of u or 2J is in the image of u'. Notice that if 2 J 
is in the image of u then Pf C u = 0. If 2 J is in the image of u' then u'(2J — 2M) — 2 J and 
hence u = t for some t € ^2M- Thus we may replace the sum over m (|7.5|) with a sum 
over 3\ M . Consequently, 

F N (^;s)= J2 sgn(u)PfiV-PfC u 

(L.M) uelti, 

J 

= E E sgn(u)Pfi? u , -PfC7 u , (7.6) 

M=o tte3 ^ 

where the second equation follows since the summand has been made to be independent of L. 
The final step in the proof of Theorem 11.141 will be establishing the following lemma. 

Lemma 7.7. Suppose that R and C are antisymmetric 2 Jx 2 J matrices, and let U = R+C. 
Then, 

J 

PfU=Y^ Yj sgn(u)Pfi? U ' • PfC u . 

It follows that F N {<&; s) = Pf(i? + C). From the definition of Up we see that Up = R + C, 
and hence 

F N ($;s)=PiUp. 

7.3. The Proof of Lcmma\L% 

Instead of computing the Jacobian of £?£,m we will compute the Jacobian of the map E' L M : 
R L x C M -> R N given by E' LM (a,/3) = b where 

N L M 



]T b n x N - n = \{{x + a t ) Y[(x + I3 m )(x + /3 m ). 



n=l 
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That is, E' L M (a,(3) = El,m{—ol, —ft). Clearly Jac E L , M (a, P) = 2&cE' LM (a, f3). 
The nth coordinate function of E' L M (a, (3) is given by, 

e„ = e n (a,f3) = e n (ai, . . .,a L ,0i,0i, ■ ■ . ,0m, (3m), 

where e„ is the nth elementary symmetric function. We will use the standard convention that 
e = 1 and e„ = if n < 0. We also specify that if 1 < I < L then e n j = e„ ; £(a,/3) is the 
n-the elementary symmetric function in all of our variables except ct(. Similarly if 1 < m < M 
then we define e' n m = e' n m (a 7 (3) to be the n-th elementary symmetric function in all of our 
variables except m and (3 m . 

Using these definitions it is easy to see that for 1 < I < L, 

de n 

Setting /3 m = x m + iy m we may compute the partial derivatives of e„ with respect to x m 
and y m . We may categorize the monomials in e n into four types: those which contain /3 m but 
not (3 m , those which contain (i m but not f3 m , those which contain both (3 m and (3 m , and those 
which contain neither (3 m nor (3 m . That is, 

e n = {Pm + 0m) e n -l,m 0m0m e n-2,m + e n,m- 

Or what amounts to the same thing, 

e n = 2x m e n _ l m + (x m + y rn ) e n _ 2 m + e n m . 

It follows that 



de n 



and 



de r , 



3x m n - L ' m "^n-Z,m Q y 

The Jacobian of E L}M is thus | det J| where the jth row of J is given by 

( 2 4-i,i + 2x i 4-2,1 2 yi e j-2,i ••• 2e 'j-i,M + 2 ^m 4_ 2;M 

e^— 1,1 e j-i,2 

Now let / be the L x L identity matrix and let 

B=U 1 1 

Define A to be the N x N block diagonal matrix 

(B 

B 

A 



2y M e' 



■j-2,M 

ej-i,L)- 



B 



V 



I) 



and set J' = J A. The j-th row of J' is given by 

(4-i,i + & 4-2,1 + 01 4-2,i 

4-1. M @ M e 'j-2-M 

e 3-l,l 

and it is easily seen that | det(J)| = 2 M | det(J')|. 

Now, let 1 < I < L and define fe(x) to be the polynomial 

L M 

Mx) = n ( x + n> + ' 3 ™)( x +m= 



4-1, M + /^M 4-2 



(7.7) 



e J-l,2 



s j-l,U 



fc=i 



AT 
n=l 
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Similarly, for 1 < to < M define g m and g m by 



M 



9m(x) = Y[(x + a e )l (x + f3 m ) Yl (x + f3 k )(x + (3 k ) 

N 



n=l 



and 



M 



9m(x) = ]J(x + a e )l (x + f3 m ) Y[ (x + p k )(x + ft) 

N 

= J^( e n-l,m + fi-m, e n-1,m) x 



n=l 



Notice that the coefficient vectors of the g mi g m and fe appear as the columns of J'. This is 
useful in light of the following orthogonality relations. By construction, fi(—j3 m ) = fe(—f3 m ) = 
for all 1< £ < L and 1 < to < M, and 



M 



M-a k ) = 



= J \{(- a t + a J )\{{-a l + (3 m ){-a t + [3 m ) if k = I, 

otherwise. 







m— 1 



Similarly, g m (-a e ) = g m (-ae) = g m (-Pm) = 9m(-Pm) = for all 1 < I < L and 1 < to < M , 
and 

1=1 V j^m ) 





5m (-ft) 

and 



TO, 



otherwise, 



Y[(-/3 m + a e ) (-/3 m + ft„) J] (-/?m + ft) (-/3 m + ft) if fc = to, 

£=1 V i#m / 





g m (-Pk) = < 

Now, let D be the N x N matrix given by 



otherwise. 



D = 



n- 


-Wi) N ~ 


l 


(- 


-Wi) N - 


-2 


"ft 


1\ 


(- 


-pi) N - 


l 


(• 


-pl) N - 


-2 


• "ft 


1 


(- 




-l 


(- 




-2 


• ~ft7 


1 


(- 




-l 


(- 


-(3m) N 


-2 


—/3m 


1 


(- 


- ai ) N - 


-l 


(" 




-2 




1 


(- 




-l 


(" 


-a 2 ) N ~ 


-2 


— «2 


1 


\(- 




-l 


(- 




-2 




V 



Clearly, Z? is a permutation matrix times the N x N Vandermonde matrix in the variables 
-ft, -ft, . . . , -/3m, ~/3m, — ai, — a2) • ■ • , —oil- 
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And thus, 

\detD\ = \detV a ' \ . 
Now, from the definitions of D and J' (equation 17. 7|) we find 



DJ' = 



(gi{- 


-pi) 


9i(- 


-Pi) 


9m{- 


-Pi) 


9m{- 


-Pi) 


h(- 


-Pi) 




-Pi)\ 


9i{- 


~Pi) 


9i(- 


-Pi) ■ 


■ 9m{- 


-Pi) 


9m(- 


-Pi) 


M- 


-Pi) ■ 


■ h{- 


-Pi) 


9i(- 


-M 


9i(- 


-Pm) ■ 


■ 9m(- 


-Jm) 


9m{~ 


-~Pm~) 


h(- 


-Wi) ■ 


■ h(- 


-Pm) 


9i(- 


-Pm) 


9i(~ 


-Pm) 


9m{- 


-Pm) 


9m{~ 


-Pm) 


fi(- 


-Pm) 


h(- 


-Pm) 


9x{- 


-Oil) 


9i(- 


-Oil) 


9m{- 


-Oil) 


9m{- 


-ai) 


h(- 


-ai) 


h{- 


-ai) 


9i(- 


-OL2) 


~9i(- 


-a 2 ) ■ 


■ 9m(- 


-a 2 ) 


9m{- 


-a 2 ) 


h(- 


-a 2 ) ■ 


■ h{- 


-a 2 ) 


\9i(- 


-oll) 


h(- 


-oil) ■ 


■ 9m(- 


-a L ) 


9m{- 


-a L ) 


h(- 


-oil) ■ 


■ h(- 


'-oll)) 



But from the orthogonality relations we see that this is in fact a diagonal matrix, and 

L M 

llM-at) II 9i(-P~i)~9i(-Pi) 



| det(DJ') 



rn—1 



= \detV a ' p \ 



But this implies that | det J'\ = \detV a ' ,3 \, and hence 

Jac(-E i)Af (c*,/3)) = | det J| = 2 M |det./'| = 2 M \detV a '^\ 

7.4. The Proof of Lemma\7^ 

Applying (|6.5p to 7 = (cc, (3), we see that 



j<k 



L M 

-MI II 

1=1 m=l 



ad 



M 



J] \Pn - Pm\ 2 \Pn-Pm\ >]] 2&(P m ). 



And hence, 



M 



\detV a - p \ = {-i) M \ n sgn(a fc - a 3 ) J] sgn3(/3 m ) ) detV a ' p 

j<k m—1 



(7.8) 



As in the end of Section ^. II we may replace the monomials in the Vandermonde matrix with any 
complete family of monic polynomials without changing its determinant. That is, detV a, P = 
det W a, P. Using the Laplace expansion of the determinant (|7.1j) with u = ie ^2Mi we see that 

detW aJ3 = J2 Bgn(t)drtW£'"-detW££. 

ltJ 2M 

Notice that the minors of the form consists of elements from the first 2M columns of 

W a 'P. These columns are not dependent on a. and thus we may write these minors as Wf t . 
Similarly we may write the minors of the form W"f as W? v . It follows that 

det V afi = s & n ( l ) det W i,i ■ det W v,v » ( 7 - 9 ) 

and the Lemma follows by substituting (|7.9p into (|7.8[) and simplifying. 
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7.5. The Proof of Lemma El 
We start by setting 

C l 



1 

L! 



detW { ? t , -PfT a m c/>( ai )- s \ d\ L {a) 



(7.10) 



e=i 



where t is an element of H^m- Expanding det W" t , as a sum over Sl allows us to write O as 



<jes l m U=i ) U=i 



A(fc)K W ) } WT a d\ L (a). (7.11) 



Recalling that for each a £ Sl, PfT a ' a = sgn(cr) Pf T a , we use the change of variables a 
o~ x ■ ol to write © as 



Sgn(er X ) 



f[^M~ fl i Il P t(fc)K) PfT«dA L (a) 



,fe=i 



Substituting this into (|7. 1 1[) we see that the sum over Sl exactly cancels 1/L\. That is, 

|n^)" Sp tw(^)| PiT a d\ L (a). (7.12) 
Using Lemma I7T51 and setting K to the integer part of (L + l)/2, we may write Pf T a as, 
PfT c 



XI SgI1 ( T ) ) II SgI1 ( a ^(2fe) - «r(2fc-l)) > 

' ren 2K U=i J 



Substituting this into (|7.12[) we find 

( L 



ren 2 



K 



° = ^7 Z~2 s S n ( T ) , "I II Sp tw( a ^) f -| II s g n ( a r(2fc) - ar(2fe-l)) f dA L (a), 



,fc=l 



© 



(7.13) 

If L is odd, then for each r G LT2if there is a fc such that a T (2fc ) = oll+i- If we set l — 
r(2ko — 1) then we may write © as 

© = 0(a^o)~ Sp t'(^)(a£„)| Yl 4>{ a e)~ Sp t'(£)( a e)^j^ Jl 8 S n ( a r(2k) - av(2fc-i))j 

^ fc=i 

x^ > (t'oT)(2fe)(Q<r(2fe))-P(t'oT)(2fe-i)(aT(2fc-i)) sgn(a r ( 2 fc) - a T (2fc-i)) j, 

where the second equation comes from reindexing the first product by £ i— > r _1 (£) together with 
the fact that 2(if — 1) = L — 1. Substituting this into (|7.13p and applying Fubini's Theorem 
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we find 



° = JO_ E s S n ( r ) ^ Sp (i'°r)(2k -i){x)dx 
' ren 2K R 



K 

n 



(x) s (j)(y) s P(vo T )(2k)(y)P(t' r)(2k-i){x)sgn(y - x)dxdy 



4^ ^(r)\f[(P { 

' T6n 2Jf I *=i 



'(t'o T )(2fe-l)! P(t'oT)(2k))l 



(x) S P(i' 0T )(2k ~l)(x)dx. 



The use of Fubini's Theorem is justified since the integrals in the latter expression converge if 
3?(s) > N. Recalling the definition of t' gives (t' Q o r)(2k ) = 2 J, and hence 

1 K 
° = JF\ E s § n ( r ) R ^ ~ !)» T ( 2fc o)] II ^ K 2fc - 1). T ( 2fc )l ( 7 - 14 ) 



^=1 



Similarly, when L is even, O is given by 

(x)~ s Hvy S P{t'or)(2k)(y)P(t'oT)(2k-i)(x) sgn(y - x) dxdy 

K l l \ 1L \* yvuT)(2k-l),P(t>OT)(2k))»\- ( 7 -15) 

Regardless if L is even or odd, (|7.15p and (|7.14|) imply that, 

1 £ 
°=ln E sgn(r)[]^K2fc-l),T(2fc)]=Pfi? t ,. 



^ e ■**(*•){ n 

' ren 2K I fc=i ' 

h e ^(r){n^ 

' ren 2K I fc=i 



ren 2/ - 



fc=i 



7.6. The Proof of Lemma\7M 
To prove Lemma 17.61 we set 



H) 



M! 



From the definition of Wf, we can write 



A I 



detTy^ t = ^2 S g n ( T ) [I P (t«T)(2m-l)(/3m)i'(tor)(2ra)(/'m) 
reS2M \,m=l 



Substituting this into we see 







t£S 2 h ^ l.m=l 



;)-^(/3 m )- s sgnQ(/3 m ) 



A/ 



X |JJ P(to T )(2n-l)(/9n)-P(ior)(2n)(/9n) j d\ 2 M((3). 
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When 3?(s) > N this integral converges, and hence we may use Fubini's Theorem to write 

M 



2 M Ml 

tES^m K m=l 



h*\ e sgn ^i n(- w ) 

X f(tor)(2m-l) (/3)-P(toT)(2m) (/?) Sgn 3(/3)cZA 2 (/3) 
E S § n ( T ) II ( P (tor)( 2m -l),^(tor)(2m))c, 

tES 2 m m=l 

which is Pf C t . But, by definition, t = t , and hence = Pf C to as desired. 



7.7. The Proof of Lemma [771 

Before proving Lemma 17.71 we present two alternative formulations of the Pfaffian. Indeed, 
many authors give one of these two formulations as the definition of the Pfaffian. 
An alternative proof of Lemma 17.71 is given in [11] . 

Lemma 7.8. Let U be a 2 J x 23 antisymmetric matrix. 

(i) Let n 2 j denote the subset of S2 j composed of those o with a(2j) > a(2j — 1). Then, 

1 J 
WU=- sgn(r)l[U[T(2j-l),T(2j)}. 

' Ten 2J 3=\ 

(ii) Let Vi,V2, ••• , v 2 j be the standard basis for R 2J , and let uj be the 2-form given by 

uj = Y U \j> fc ] v i A v fc- Then, 

j<k 

uj A uj ■ ■ ■ A a; = Pf C7 ■ Vi A v 2 A • • • A v 2 j. 

,7 

Proof. Let G 2 j be the subgroup of S^j generated by the transpositions (2j — 1 2j) for 
7 = 1,2,..., J. Then if 7r G G 2 j and r € S^j, the antisymmetry of U implies that 

J J 
sgn(r) Y[ U[t(2j - 1), t(2j)} = sgn(r o tt) C/[(r o 7r)(2j - 1), (r o 7r)(2j)]. 
3=1 3=1 

It follows that we may replace the sum over S2J m the definition of the Pfaffian with a sum 
over left cosets of G 2 j. Each coset contains 2 J elements and II 2 j forms a complete set of coset 



representatives which establishes (i) 
we write 



To prove (ii) 



lc i A^A__A 5 =l| E •■■ E (n^' m ^) • A( v ^ Av ^)|- 

j K fci<mi fc 2 <iTi2 kj<mj \i=l / 1=1 ) 

Notice that if any two of the indices of summation are equal, then the summand on the right 
hand side of this expression is identically zero. Thus we may replace the J-fold sum with a 
single sum over LI 2 j to write 

-L yA^A---A^ = ji [ Y[U[T(2j-l),T(2j)]) ■ A( v -(2«-i) Av ^))|- 
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(The big wedge notation is unambiguous here since this wedge product is independent of order). 
Then, (ii) follows from (i) by noting that 

J 

f\ (v T ( 2 £-i) A \ T (2t)) = sgn(r) ■ \ x A v 2 A • • • A v 2J . 



Proof of Lemma 1 7. 71 Following (ii) in Lemma 17.81 define the 2- forms g and x by 
g = m]v fc A v m and x = 2J C[k, m]v fe A v m , 



where v l7 v 2 , . . . , v 2 j is the standard basis for M. 2N . Then 
j , j 



X A x A • • • A x A g /\ g f\ ■ ■ ■ f\ g 
v „ ' v , ' 

j J-j 



j=0 ^ J 



XAxA---Ax M 



(J-jY- 



g A p A • • • A g 

" v ' 

J-j 



(7.16) 



Using the Lemma 17.81 and the linearity and alternating property of the wedge product it can 
be established that 



- x A x A ■ ■ ■ A x = Pf Cu ■ v "(!) A v "(! 



) A • ■ ■ A v u(2j) 



An analogous formula holds for 1 / {J — j) \ • g f\ • • • f\ g. Substituting these expressions into (|7.16[) 
we find 

1 J 

j]/\U + $ = Yl ll Pf C„ ■ Pf i? t ■ v u(1) A ■ ■ ■ A v u(2j -) Av l( i)-Av l(2J _ 2)) . (7.17) 



j'=i 



It is immediately clear that the summand is identically zero unless the ranges of u and t are 
disjoint — that is, unless t = u'. Thus we may remove the sum over 3 2 j_ 



j { ACx + e) = i Yl PfC u -PfiJtt'-v 4u(1) av 1u 



2,; ■ 



(2) A • ■ ■ A V, u ( 2 j) 



where we combined the two wedge products in (|7.17p into a single product using the definition 
of l u . Now, since sgn(u) = sgn(t u ), 

ji/\(x + Q)= \ J2 sgn(u)PfC u -Pfi? u < 1 - Vl Av 2 A---Av 2J . 



And the lemma now follows from Lemma 17.81 
7.8. The Proof of Corollary\TjM 



□ 



Suppose that P,Q £ R[x] are either both even or both odd. Then, if (j>(—(3) = </>(/3) for every 
£ C it is easy to verify that (P, Q)m = (P, Q)c = 0. Notice if P £ R[i] is odd, then 

4>{x)~ s P(x) dx = 0. 

Corollarv ll.lGl is a consequence of the following lemma. 
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Lemma 7.9. Suppose that U is a 23 x 2J antisymmetric matrix such that U[j, k] = if 
(j - k) ee mod 2. Then, 

VfU = det A, 

where A is the J x J matrix whose j, k entry is given by A[j, k] — U[2j — 1, 2k]. 
Proof. 

1 J 

PfC/ = Fj! E sgn(a)l[U[a(2j-l),a(2j)}. (7.18) 

Clearly the product in this expression is except for permutations a such that 

a(2j - 1) - a(2j) ee 1 mod 2 for j = 1, 2, . . . , J. (7.19) 
Let Gat denote the subgroup of Sn given by 

Gjv = {cr e S'at : (o-(n) - n) = mod 2, n = 1, 2, . . . , N}. 
Given a e Sn satisfying (|7.19p . define cr* e Gjv by 

<r(2j) if a(2j) is even, 



**(2j) 
and 



er(2j - 1) if cr(2j) is odd, 



f (7(2 j) if a(2j) is odd, 
1 J J ~ \ a(2j - 1) if a(2j) is even. 



Notice that a and a* differ only by a product of transpositions of the form (2j — 1, 2j) where 
7 = 1,2,..., J. Moreover, since U is an antisymmetric matrix, 
J ./ 
S gn(<7) U[a(2j - 1), ( t(2j)] = sgn(a*) [] C/[<7*(2j - l),a*(2j)]. 

Clearly, each a* 6 Gat represents 2 J different permutations cr £ Sat each of which satisfy 
(|7.19p . We may thus replace the sum over Sn in (|7.18|) with a sum over Gn to find 

Pf ^) = 7! E Bgn(a*)]JU[a*(2j-l),a*(2j)}- 
' <r*eGjv j=i 

Now, since elements oi Gn permute even integers and odd integers disjointly we have Gn is 
isomorphic to Sj x Sj, and hence 

1 - 

tBSj ttSSj j'=1 

But, by [9l Lemma 3.1] this is exactly det A. □ 
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